Computation of geometric measure of entanglement for pure multiqubit states 
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We provide methods for computing the geometric measure of entanglement for two families of 
pure states with both experimental and theoretical interests: symmetric multiqubit states with 
non-negative amplitudes in the Dicke basis and symmetric three-qubit states. In addition, we study 
the geometric measure of pure three-qubit states systematically in virtue of a canonical form of 
their two-qubit reduced states, and derive analytical formulae for a three-parameter family of three- 
qubit states. Based on this result, we further show that the W state is the maximally entangled 
three-qubit state with respect to the geometric measure. 

PACS numbers: 03.65.Ud, 03.67.Mn, 03.67.-a 



I. INTRODUCTION 

Quantum entanglement, which was first noted by Ein- 
stein and Schrodinger Q], 0] , has been extensively studied 
in the past 20 years [3j. In particular, multipartite en- 
tanglement has attracted increasing attention due to its 
intriguing properties and potential applications in quan- 
tum information processing. 

The importance of multipartite entanglement can be 
illustrated in two aspects. In respect of application, 
graph states, prominent examples of entangled multi- 
qubit states, are a useful resource for one-way quantum 
computation 4\ and fault-tolerant topological quantum 
computation [5[. Multipartite entangled states, such as 
GHZ states, are essential resources for quantum secret 
sharing @, Q. In addition, multipartite entangled states 
can serve as multiparty quantum channels in virtue of 
teleportation Q. In respect of theoretical interests, mul- 
tipartite states display stronger nonlocality, one of the 
key features of quantum physics 0-[H| • Quantum cryp- 
tography beyond pure entanglement distillation has been 
generalized to multipartite bound entangled states [l^ . 
What's more, recent progress in experiments has made 
accessible more multipartite entangled states, such as the 
GHZ states [Hj|, W states [l4[, six-photon Dicke states 
[l5l - fl7| etc. Methods for detecting such states have also 
been developed [l8j . 

Given an entangled state, a natural question to ask 
is how much entanglement is contained in this state. 
In quantum information theory, entanglement is usually 
quantified by entanglement measures [19J. An entangle- 
ment measure is an entanglement monotone, which can- 
not increase under local operations and classical com- 
munications (LOCC), and equal to zero for only clas- 
sically correlated (separable) states |20|. Hitherto, the 
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most well-known entanglement measures arc defined for 
bipartite states, such as entanglement cost and distillable 
entanglement [2fJ,[21|. For pure bipartite states, there is 
essentially a unique entanglement measure, the von Neu- 
mann entropy of each reduced density matrix, which is 
easily computable (22| . 

For multipartite states, while a lot of entanglement 
measures have been proposed (3l. [23l - |25j . the characteriza- 
tion of multipartite entanglement is far from being com- 
plete. It is generally difficult to calculate such measures 
even numerically. Moreover, the existence of many types 
of inequivalent entanglement defies a unique definition. 
Different entanglement measures often induce different 
orders and even lead to different maximally entangled 
states. For example, the Bell state |*) = ^=(|00) + |11)) 
is the maximally entangled state of a two-qubit system 
for all measures, since it violates the Bell inequality most 
strongly. However, its multipartite analog, the GHZ 
state | GHZ) = ^(|000) + | 111) ), is maximally entangled 
only under some specific entanglement measures, such as 
three-tangle [26|, [27| ■ It is not a maximally entangled 
state under the definition in Ref. (28|, and the geomet- 
ric measure of entanglement, which is one focus of the 
present paper. 

On the other hand, some geometrically motivated mul- 
tipartite entanglement measures have been providing us 
insights on quantum entanglement. One prominent ex- 
ample is the geometric measure of entanglement (GM), 
[23l l2~ij which quantifies the minimum distance between 
a given state and the set of product states. In addition 
to providing a simple geometric picture, GM has signif- 
icant operational meanings. It is closely related to opti- 
mal entanglement witnesses 0, [2!| , and has been shown 
to quantify the difficulty of multipartite state discrimi- 
nation under LOCC [30]. Recently, GM has also been 
applied to show that most entangled states are too en- 
tangled to be useful as computational resources [3l|. In 
condensed matter physics, GM has been utilized to study 
the ground state properties and to characterize quantum 
phase transitions 32- 3o| . 
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There have been extensive literatures on the quantita- 
tive calculation of GM for both pure and mixed states 
[HH^ISSl- The qualitative analysis on GM has also 
received much attention [U |42j ■ In addition, a few nu- 
merical methods have been developed for computing the 
GM of multipartite states, such as the algorithms pre- 
sented in Refs. [43|, [44|, which allow repeated analyt- 
ical maximization according to a subset of the param- 
eters with a high efficiency. However, our knowledge 
about GM is still quite limited. Even for pure three- 
qubit states, there is no complete analytical solution. In 
addition, it is still uncertain which state is the maximally 
entangled with respect to GM, although the authors of 
[27l ] conjectured that the W state is such a candidate. 
Thus it is desirable to compute GM analytically for more 
entangled states, which is another focus of the present 
paper. 

In this paper, we would like to compute the GM for 
several families of multipartite pure states and determine 
the maximally entangled three-qubit states with respect 
to GM. Throughout the paper, by symmetric states, we 
mean those states which are supported on the symmetric 
subspace of the whole Hilbert space. 

First, we present an analytical method for comput- 
ing the GM of symmetric multiqubit states composed of 
Dicke states with non-negative amplitudes by virtue of 
a recent simplification on GM of symmetric states [42j . 
Next, we analytically compute the GM for symmetric 
three-qubit states. Combining with the results in Ref. 
[36j | , we provide a complete analytical solution to GM of 
any symmetric pure three-qubit states. Recall that many 
important multiqubit states accessible to experiments so 
far are symmetric, e.g. GHZ states W states flij . 



and Dicke states [lg, [l7[ etc. Our results may hopefully 
help analyze these states in experiments. 

Second, we introduce a canonical form of pure three- 
qubit states based on the canonical form of two-qubit 
rank- two states developed in Ref. (45[. In virtue of this 
canonical form, we study the GM of pure three-qubit 
states systematically and derive explicit analytical for- 
mulae of GM for a three-parameter family of three-qubit 
states. Starting from these results, we prove that, up to 
local unitary transformations, the W state is the unique 
maximally entangled pure three-qubit state with respect 
to GM, confirming the conjecture made in Ref. [27| . 

The rest of the paper is organized as follows. In Sec. 

II, we propose analytical methods for computing the GM 
of symmetric multiqubit states with non-negative ampli- 
tudes and that of symmetric three-qubit states. In Sec. 

III, we derive analytical formulae of GM for a three- 
parameter family of pure three-qubit states and prove 
that the W state is the maximally entangled state under 
GM. We conclude in Sec. IV. 



II. ANALYTICAL METHOD FOR COMPUTING 
GEOMETRIC MEASURE (I): SYMMETRIC 
STATES 

The definition of GM of bipartite or multipartite pure 
states is motivated by the following simple geometric idea 
of entanglement quantification: the farther away from the 
set of separable states, the more entangled a state is [23| . 
Given a pure state \ip) of a joint system composed of 
subsystems A, B,C, ■ ■ ■ , define G(\ip)) as the maximum 
overlap between and the set of product states, that 



G(m := max |(^)|, 

\tp) = \a)\b)\c)... 



(1) 



where the normalized one-particle states \a), \b), |c), • • ■ 
belong to subsystems A, B,C,. . ., respectively. G(\tp)) is 
manifestly invariant under local unitary transformations. 
It obtains the maximum value 1 only for product states 
and is thus an inverted entanglement measure. The GM 
of a pure state is defined as follows: 



£ G (|^»:=l-G(h/>» 2 , 



(2) 



or in another version —2 log G(\ip)) sometimes. In this 
paper we will follow the definition in Eq. @. It can 
be extended to the GM of mixed states by convex roof 
construction [23| according to the same idea as in the 
definition of entanglement of formation (20| : 



K 



G (p) := min V" 

p=]CiP«IV>i>W'il . 



(3) 



Eq(p) has been shown to be an entanglement monotone 



by T.-C. Wei and P. M. Goldbart |23(. An alternative 
definition of GM for mixed states will be introduced in 
Sec. Mil in a different context. 

Clearly, E G {\^)) in Eq. © is determined by G(\ip}) 
in Eq. ([T]). From now on we focus on G(|t/>)) of pure 
states and call it GM too, if there is no confusion. 
For a pure bipartite state, the GM is equal to its largest 
Schmidt coefficient. The problem becomes difficult for 
pure multipartite states, since there is no Schmidt de- 
composition in general. The difficulty lies in the linearly 
increasing number of optimization variables parametriz- 
ing the product states \a) \b) \c) ■ ■ ■ in Eq. (fT]), as the num- 
ber of parties increases. In fact, only a few partial results 
are available on this problem ]27l [36ll37j. 

Recently, the authors of Ref. [42[ proved that, for a 
symmetric pure state |'0 sy ), it suffices to consider sym- 
metric product states in the maximization in Eq. ([TJ, 
that is, 



G(|V sy » 



max 

= |o>|a>|o> 



(4) 



This result can greatly simplify the calculation of GM for 
symmetric states. In the rest of this section, we derive an- 
alytical solutions for two families of states, respectively, 
in virtue of this result. In Sec. lII Al we analytically derive 
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GM for symmetric multiqubit states with non-negative 
amplitudes in the Dicke basis. In Sec. Ill B[ we derive the 
analytical solution of GM for symmetric pure three-qubit 
states based on the previous work (36|, thus solving this 
problem completely. 



Unfortunately, the present method does not apply to 
arbitrary symmetric multiqubit states, e.g., those states 
in Eq. ([5]) having negative or complex amplitudes. For 
more complicated states, numerical methods are indis- 
pensable for computing their entanglement measures, 
see, for example, Refs. [43|, |44|. 



symmetric multiqubit states with non-negative 
amplitudes 



B. symmetric three-qubit states 



In this subsection, we compute the GM for pure sym- 
metric multiqubit states with non-negative amplitudes 
in the Dicke basis. More explicitly, we investigate the 
A^-qubit state 



|^ symq ) 



N 

E< 

rn— 



\m,N), 



(5) 



where a m > 0, and \m,N} is the Dicke state [15| defined 

as 



\m,N) := 




(6) 



where P^s denote the set of all permutations of the spins. 
By definition, Dicke states are symmetric; so the state 
|'0 symq ) is also symmetric, and we can apply Eq. (|4|) to 
computing its GM. Let \a) — cosa|0) + e lS sina|l) with 
a £ [0, §] and 9 € [0,2tt]; then the GM of the state in 
Eq. (jSJ) reads 



G(|V symq }) 



max 

|^> = |a>|o)|a>. 



< max y 

\<p)=\a)\a)\a,)...£± Q 



max > 



m=0 x 

N ^12 



a m cos 



N- 



a m cos N m a sin m a 



cim cos N m a sin m a, 



(7) 



where the equality holds when 9 = 0. Equation ([7]) con- 
tains only one variable a, so one can easily find out the 
maximum. For example, let x = tana, then one can 
convert G 2 (|'0 symq )) into a rational fraction A(x)/B(x), 
where A(x) and B(x) are both polynomials on x. By 
calculating its derivative we can find out the maximum 
in Eq. explicitly. 

A similar idea can be applied to calculating the GM of 
any symmetric multi-qudit state with nonnegative am- 
plitudes in the generalized Dicke basis; again the number 
of free variables can be reduced by half. Recently, the 
additivity of GM of states with non-negative am plit udes 
was proved, i.e., G(\a) ® = G(\a)) G(|/3)) [U. So 
we can compute the GM of |?/' symq ) <Ei \(3) if we know the 
GM of IS) too. 



In this subsection, we compute the GM of symmetric 
pure three-qubit states. Such states can always be con- 
verted into the following form with suitable local unitary 
operations [36j : 

|$) = ,g|000) + t(|011) + 1 101) + 1 110)) + e^h\ni), (8) 

where g,t,h > and 7 € [—§,§]■ So it suffices to calcu- 
late the GM for the state |$). 

Analytical formula of the GM is already known if at 
least one of the three par ameters g, i, h is vanishing, or 

7 = o,±§ m, m ni,i33, 

the family of states with 



Hence, we can focus on 



g,t,h>0, g 2 + 3t 2 + h 2 = l, T € (--,0)U(0,-). 

(9) 

The authors of Ref. [37j have reduced the task of com- 
puting the GM of the state |$) to solving the following 
system of equations of the three variables <p, 9, A (see also 
appendix A): 

2ht cos 7 + 2t(g + t) sin 9 cos ip~ 

2ht cos 7 cos 9 =A sin 9 cos <p, 
(10) 

2ht sin 7 — 2t{g — t) sin 9 sin (p— 

2ht sin 7 cos 9 =A sin 9 sin <p, 
(11) 



(g 2 -t 2 )(l + cosf 



K 2 (l 



2ht cos 7 sin 9 cos ip — 2ht sin 7 sin 9 sin ip =A cos 9. (12) 

For each root (ifj, 9j) of Eqs. ((T0)) - p2)) . we can obtain a 
GM candidate of the state |<&) via the following formula, 
according to Eq. (|32|) in appendix A, 

G 2 (|$)) = i [3 - 2t 2 + 4(1 - 2h 2 - At 2 ) cos 9 j 

+ (1 — 6t 2 ) cos 29 j + igt cos 2<pj sin 2 0j 
9 ■ 9 ■ 

+ 32^cos( 7 -( < 3 J )cos^-sin 3 -|-]; (13) 



the GM is the maximum over all the GM candidates: 



G 2 (|$)) =maxG 2 (|$)). 



(14) 



We shall solve Eqs. (ITUl) - (|12p in two cases separately; 
the second case consists of three subcases. In each case 
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we obtain one or a few GM candidates by computing 
Eq. ([13")) with the solutions to the system of equations. 

Case 1. Suppose 9 = 0; then the phase ip does not 
play any role, and Eqs. (fTUf and become identities, 
while Eq. ([12")) determines A. In this case we get a GM 
candidate via Eq. (|13l) as follows: 



G\m)=g 2 



(15) 



Case 2 To satisfy Eqs. tpXi ]) — (fl^ ]> . the roots <p 
for k = 0, 1, 2, 3 lead to # = 0, which is already discussed. 
Moreover, 9 = ir cannot be a legal solution of Eqs. (fTQ|) 
or (fTT|) . so this choice is excluded. Hence, it remains to 
solve Eqs. (fTU)) -(fT2" )) under the assumption that 

^G(0,|)U(|, 7 r)U(7r,3|)U(3| I 27r), 9e (0,tt). 

(16) 

From Eqs. (fTUl - [T2")) . we can determine A as a function 
of 9 and ip. Inserting this solution into Eq. (jlpp and 
Eq. (1111) . we can obtain two equations about 9 and ip: 
eql(<p,8) = and eq2(p,9) = 0, respectively. This fur- 
ther implies that either 



tan - = — — -. 

2 fi esc zip sm( r y — pj 



tan - = — — r. 

2 h esc 2(p sm(7 + p) 



Combining either of them and eql(<y2, 9) = can lead to 
a set of solutions. 



Case 2.1 There is a simple solution tan ip 



_ t+g 

t-g 



tan 7. 



The variable 9 can be determined via either Eq. (fl7|) 
or (fTg|) . which lead to an identical result in this case. 
Inserting this solution into Eq. (fl"3|) . we get another GM 
candidate: 



Glm)=g 2 - 



(g 2 -t 2 ) 



2\3 



t 2 - 2t 4 + g 2 - 6g 2 t 2 - 2gth 2 cos 27 ' 

(19) 



Case 2.2 By combining Eq. (JTTJ) and eql(tp, 9) = 0, 
we can get two polynomial equations: 



^2c u (g,t,h,j) cos 1 2<p = 0, 

i=0 
4 

c 2l (5, t, h, 7) cos 1 2<p = 0, 



(20) 



i=0 



as well as Case 2.1, which has already been handled. 
Since Eqs. (|2"0")) are quartic equations on cos2(p, we can 
analytically derive their roots. We may obtain up to 16 
different phases ip e [0,27r]. The variable 9 can then be 
determined via Eq. (|17p. Hence, we can derive up to 16 
GM candidates Gj(|$)) for j = 3, . . . , 18 via Eq. (1131) . 



This differs a bit from Case 1 and Case 2.1, where only 
one GM candidate is given respectively. 

Case 2.3 Similar to Case 2.2, by combining Eq. (fT8|) 
and eql(</3, 9) = 0, we can get two quartic polyno- 
mial equations on cos 2ip, which are analytically solvable. 
Again we may get up to 16 GM candidates Gj(\&)) for 
j = 19, . . . , 34. This finishes the discussion of Case 2. 

Now we have all the GM candidates Gj(\Q)) for j = 
1, . . . , 34. The maximum of them is exactly the GM of 
the state |$) in Eq. ©. 

For the convenience of the readers, here we repeat the 
main steps for deriving the GM of the symmetric three- 
qubit state |<E>). 

Step 1. Compute Gi(|$)) and G 2 (|$)) via Eqs. (JT5) 
and (Ti~9")) . respectively; two GM candidates can be ob- 
tained. 

Step 2. Compute Gj(|$}) for j = 3, . . . , 18 via 
Eq. with roots (<Pj,6j) of Eqs. JT7J and (HU; up 
to 16 GM candidates can be obtained. 

Step 3. Similar to Step 2, with the roots (ipj,9j) for 



j = 19, . . . , 34 of Eq. ([18)) and quartic equations similar 
to Eq. (|2T)]) . up to 16 GM candidates can be obtained via 
Eq. (USD. 

Step 4. The maximum of all 34 GM candidates is 

( 17 ) exactly the GM of the state |$). 

In conclusion, we have provided a method for analyti- 
cally deriving the GM of the symmetric three-qubit states 
in Eq. © with 7 G (-§ , 0) U (0, |). The special cases 

(18) 7 = 0, ±77 have been addressed in Ref. (36[. Calculation 
shows that our result approaches their result when 7 ap- 
proaches these special values. Hence, we can now com- 
pute the GM of any symmetric pure three-qubit states. 



III. ANALYTICAL METHOD FOR 
COMPUTING GEOMETRIC MEASURE (II): 
MAXIMAL ENTANGLED STATES AMONG 
PURE THREE-QUBIT STATES 

In this section, we introduce a canonical form of pure 
three-qubit states based on the canonical form of two- 
qubit rank- two states developed in Ref. [45[ . By virtue of 
this canonical form, the GM of pure three-qubit states is 
studied systematically. In particular, we derive analytical 
formulae of GM for the family of pure three-qubit states 
one of whose rank-two two-qubit reduced states is the 
convex combination of the maximally entangled state and 
its orthogonal pure state within the rank-two subspace. 
Based on these results, we prove that the W state is the 
maximally entangled three-qubit state with respect to 
GM, confirming the conjecture in Ref. [271 ] . 

Our approach builds on Theorem 1 in Ref. [48| , which 
states that the GM of an n-partitc pure state \ip) is de- 
termined by any of its (n — l)-partitc reduced states p, 
that is, 

G 2 (m=g(p). 
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Here g(p) is an alternative definition of geometric mea- 
sure and has nothing to do with the parameter g intro- 
duced in Eq. ©: 

g(p)= max tr[p(pi <g> • • • <g> p n -i)] , (21) 

P\,...,p n -1 

where pi, ... , p n —\ are pure single-particle states, namely 
Pi = | ■ In addition, to any closest product state 

p\ ® ■ ■ ■ ® Pn-i °f P, there corresponds a unique closest 
product state of \ip) with p\ ® • • • ® p n -\ as a reduced 
state. A closest product state of p is any pure product 
state pi (8 • • • (8 Pn-i that maximizes Eq. (|2Tj) . From the 
definition, g[p) is a convex function of p; this property 
will be frequently resorted to later. 

Note that for a mixed state p, g(p) is not the stan- 
dard definition of the GM of p (see the first paragraph of 
Sec. HI}. Nevertheless, this alternative definition is useful 
for computing the GM of any purification of p [HI]. It 
has also many applications of its own, such as construct- 
ing optimal entanglement witnesses [23l . |29| and quanti- 
fy™^ the difficulty of state discrimination under LOCC 



A. Canonical form of two-qubit rank-two states 

In this section we introduce a canonical form of pure 
three-qubit states based on the canonical form of two- 
qubit rank-two states developed in Ref. [45[ and set the 
notations useful in later discussions. 

For a pure three-qubit state, each two-qubit reduced 
state lies on a rank-two subspace of the two-qubit Hilbcrt 
space. Up to local unitary transformations, the projector 
E onto a general rank-two subspace can be specified by 
just two parameters 71,72 [HI: 

£ = i(l + ua 3 + vt 3 + Z\o\T\ + z 2 a 2 T 2 ), 
u = cos 71 cos 72 , v — sin 71 sin 72 , 
z\ = sin 71 cos 72, z 2 — cos 71 sin 72, 

with ^tt > 71 > 72 > 0, (22) 

where a% , a 2 , a 3 are the Pauli operators for the first qubit 
and Ti,T2,t 3 are that for the second qubit. Interchange 
of the two qubits leads to 71 — > | — 72, 72 — > § — 71 ■ So 
without loss of generality, we can assume 71+72 < ?, 

72 <!• 

Any rank-two state supported on £0 can be written as 
follows: 

Prk2 = ^(S + xiSi + X2E2 + ^3^3), (23) 
where £1, £2, £3 are the Pauli operators for the rank- two 



subspace [45| : 

_ 1 . 

£1 = -(sin 71CT1 + COS72T1 + sin7 2 criT3 + COS71CT3T1), 

£2 = ^(sin 72CT2 + COS71T2 + sin7icr 2 T3 + COS72CT3T2), 

£3 = 7^{v<T3 + ut 3 - z 2 ctiTi - z x a 2 r 2 + (7 3 t 3 ), (24) 

and (xi,x 2 ,x 3 ) (satisfying x\ +x 2 +x\ < 1) is the Bloch 
vector of p r ^ 2 . 

Local unitary symmetry and complex conjugation 
symmetry play an important role in determining the be- 
havior of g(p r k2) and in simplifying its calculation. Ac- 
cording to Eqs. (|2"2"|) and (|2"4")1 , simultaneous local unitary 
transformation a 3 (S> t 3 flips the sign of £ 1 , £2 , while leav- 
ing £o,£ 3 invariant, that is, 

03 ® T 3 Y, Q . 3 a 3 ® r 3 = £0,3, 

cr 3 ® t 3 £i, 2 (73 ® r 3 = -£i, 2 ; (25) 

under this transformation, the Bloch vector of p T \ 2 
changes as follows: (x\,x 2 ,x 3 ) — > (~xi,—x 2 ,x 3 ). Com- 
plex conjugation flips the sign of £2, that is 

£* = (-1)^ 2 £ J5 (26) 

where 5j t % is the Kronecker 5 function; under this trans- 
formation, the Bloch vector of p r k2 changes as follows: 
{xi,x 2 ,x 3 ) — > (xi,—x 2l x 3 ) . As a consequence of 
these symmetries, any reasonable entanglement mea- 
sure is equal for the four states p r k2 with Bloch vectors 
(±xi,±x 2 ,x 3 ), respectively. Without loss of generality, 
we can assume x\, x 2 > 0. 

If 71 =72, the simultaneous local unitary transforma- 
tion e~ ie<T3 ® T3 rotates the Bloch vector of p r k2 around the 
x 3 axis, so g(/9 r k2) is rotationally invariant about the x 3 
axis. If 72 = 0, the local unitary transformation e 1 
rotates the Bloch vector around the x\ axis, so g(/? r k2) is 
rotationally invariant about the x\ axis. 

Up to local unitary transformations, there is a one-to- 
one correspondence between pure three-qubit states and 
rank-two two-qubit states. Hence, the canonical form of 
rank-two two-qubit states provides a canonical form of 
pure three-qubit states. Moreover, due to Theorem 1 in 
Ref. [lH and the arguments given above, computing the 
GM of pure three-qubit states can be reduced to comput- 
ing <?(p r k2) of the family of canonical rank- two two-qubit 
states in Eq. (23]) with < 72 < 71 < f , 72 + 7i < f > 
< Xi,x 2 < 1 and — 1 < X3 < 1. With this background, 
we can now study the GM of pure three-qubit states sys- 
tematically. 

B. (?(p r k2) of general two-qubit rank- two states 

In this subsection we reduce the task of computing 
3(Prk2) for general two-qubit rank- two states to a maxi- 
mization problem which involves only two free variables. 
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The number of free variables is further reduced to one for 
states /0 r k2 with x-i = 0. 



Let pi 



(! + «!• <t) and p 2 



fl 



«2 ■ t) be two 

pure qubit states with Bloch vectors s% — (a, 6, c) and 
82 = (g2>&2)C 2 ), respectively. Straightforward calcula- 
tion shows that 



tr(/9 r k2Pi <£> P2) = ^(1 + axi sin 71 + 6a; 2 sin 72 
+ CCOS71 cos 72 + CX3 sin 7! sin 72 + w ■ s 2 ), (27) 



where 



w 



a(cos7 2 sin 71 — 0:3 cos 71 sin 72) 
b ( — £3 cos 72 sin 71 + cos 71 sin 72 ) 
bx2 sin 71 + ax 1 sin 72 

CX\ cos 71 + X\ cos 72 
x 2 cos 71 + cx 2 cos 72 
CX3 + X3 cos 71 cos 72 + sin 71 sin 72 



(28) 



Given p r k2 and pi, the trace in Eq. (|27|) is maximized 
when S2 is parallel to to. According to Eq. (|2"Tj) , we have 



5(Prk2) 



f(a,b,c) 



= max tr(/9 rk2 pi (8 p 2 ) 
P1P2 



— max f(a,b,c), 
4o 2 +5 2 +c a =r 



(1 + axi sin 71 + 6x2 sin 72 + c cos 71 cos 72 
+ cx 3 sin 71 sin 72 + \w\), (29) 



where \w\ denotes the Euclidian norm of w. Thus we 
have reduced the task of computing g(p r k2) to that of 
maximizing the function /(a, b, c) on the unit sphere de- 
termined by a 2 + b 2 + c 2 — 1, which involves only two 
free variables. The contours of / are in general some 
quadratic surfaces. At the maximum of / over the unit 
sphere, the contour is generally tangent to the sphere. 
This geometric picture is useful in visualizing the closest 
product states. 

Further simplification is possible for states p r k2 with 
X2 = 0. Assuming X2 — and x\ > (recall that we 
only need to consider the case < x%, X2 < 1 due to con- 
sideration on symmetry, see Sec. IIII A[) ; then f(a,b,c) 
is an even function of b according to Eq. (|29p . In ad- 
dition, f(\a\,b,c) > /(— a, 6, c) and f(a, yl — a 2 , c) is 
nondecreasing with a for a > 0. So the maximum of 
/(a, b, c) can be obtained in the parameter subspace sat- 
isfying a > 0, b = 0. Moreover, the maximum can only 
be found in this subspace if x\ > 0. Thus the calculation 
of g(prk2) can be reduced to the optimization problem 
over the single variable c. 

For rank- two subspace with 71 = 72 or 72 = 0, this 
simplification is applicable to all states p r k2, since it is 
enough to calculate g(p r k2) for states with x 2 = due to 
the symmetry discussed in Sec. IIII Al 



The W state is the maximally entangled state 
with respect to the geometric measure 



According to the discussion in Sec. IIII At to determine 
the maximally entangled states of three-qubit with re- 
spect to GM, it is enough to determine the global mini- 
mum of g{pi±2) over the set of canonical two-qubit rank- 
two states. Due to the convexity (cf. Eq. (|2"Tjl ) and 
the symmetry of <?(p r k2), given 71,72,^3, the minimum 
of g(/0rk2) (as a function of x\ and X2) is obtained at 
x\ = X2 = 0. So the global minimum of g(p r k2) can be 
obtained at states with this property. Recall that the 
state p r k2 with x\ — X2 — 0, x$ = — 1 is the maximally 
entangled state in the rank- two subspace [Hj]. Hence, 
these states are convex combination of the maximally 
entangled state and its orthogonal pure state within the 
rank-two subspace. They are interesting for a couple 
of reasons. First, the two-qubit reduced states of many 
important pure three-qubit states, such as the W state, 
GHZ state, are among this family of states. Second, from 
the result on this family of states and that on pure states, 
both an upper bound and a lower bound for (?(p r k2) of any 
state in each rank-two subspace can be obtained by virtue 
of the convexity and symmetry properties of g(/o r k2)- 

Hence, in order to find the maximally entangled three- 
qubit state with respect to GM, it suffices to investigate 
<7(/0rk2) of the states p r k2 with x\ = X2 =0. After some 
elementary algebra (see Appendix B), we derive a simple 
analytical formula of g(p r k2) of this family of states. To 
emphasize its explicit dependence on the three parame- 
ters £3,71,72, we write 3(23,71 ,72) for g(p rk 2) ■ 



gfa, 7i; 72) 



1 



(l-x 3 )[l+cos(7i+7 2 )] 
4 ' 

(1— x|) sin 71 cos 72 (cos 72 sin 71 —x 3 cos 71 sin 72) jj- 
— 2 ^| — (cos 72 sin 71—2:3 cos 71 sin 72 ) 2 J 



(30) 



( 1+Z3 )[l+cos(7i— 72)] 
4 



where I, II, III denote three intervals, I : 



II 



„(3) 



< X3 < X 



(4) 



III 



< X 3 < 1. 



III 



l<x 3 <x { 3 3 \ 
Here x$ and 



are given by 



4 M) (7i,72) := 
— sin 71 {± sin 71 



_ + [cos 71 cos 72 + (sin 71 ) 2 ] sin 72} 
1 + cos 71 {cos 72 — shi72[cos7i sin 72 + (shi7i) 2 tan 72]} ' 

(31) 

and obey the inequalities: — 1 < £3 < < x^ < 1. 

Figure [T] shows 3(^3, 71, 72) as a function of X3 for sev- 
eral different values of 71,72. 3(^3,71,72) is equal to 
I [l + cos(7i =F 72)] at £3 = ±1. This is consistent with 
the well-known result on GM of two-qubit pure states, 
which is a function of the concurrence; recall that the 
concurrence is equal to sin(7i =p 72) for the two states 
with £3 = ±1 [43. 3(^3,71,72) is equal to | at 2:3 = 0, 
which is independent of the other two parameters. This 
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/ / 







'-1.0 -0.5 0.0 0^5 1.0 



X3 

FIG. 1: p(i3,7i, 72) as a function of x$ for several different 
two-qubit rank-two subspaces. (a) 71 = 72 = 0, (b) 71 = 
§, 72 = 0, (c) 7l = f , 72 = f , and (d) 7l = 72 = f. 
Curve (a) and curve (c) coincide in a large interval, because 
71 — 72 = \ for both the rank-two subspaces, see Eq. (|30l) . 




about xz\ the minimum can be found after solving this 
equation. In particular, — 1 < X3 < at the global min- 
imum of 3(2:3, 71, 72). Figure [2] shows the dependence 
of the minimum of 5(2:3,71,72) over X3 on 71,72; the 
minimum is also the minimum of <?(p r k2) in the rank- two 
subspace. According to the figure, the global minimum 
of 3(2:3,71,72) is obtained in the rank- two subspace with 
7i = 72 = f ■ 

To determine the maximally entangled multipartite 
states is a highly nontrivial task, since it usually involves 
a massive optimization process over a large parameter 
space. Even for three qubits, the maximally entangled 
state with respect to GM is not known for sure, although 
it has been conjectured with strong evidence that the 
W state is such a candidate 27]. As an immediate ap- 
plication of the above results, we prove this conjecture 
rigorously in Appendix C. 

Theorem 1 Up to local unitary transformations, the W 
state is the unique maximally entangled pure three-qubit 
state with respect to GM. 

Theorem [T] can be generalized to mixed states according 
to the convex roof definition in Eq. 

Theorem 2 The W state is the maximally entangled 
state among all three-qubit states with respect to GM. 

IV. CONCLUSIONS 



FIG. 2: The minimum of 3(3:3,71,72) over x% in each rank-two 
subspace; the global minimum is obtained at 71 = 72 = \. 

observation implies that G 2 (\ ip)) = \ for any pure three- 
qubit state which is maximally entangled under some bi- 
partite partition; a perfect example is the GHZ state. 
Once the value of 3(2:3, 71, 72) is specified at the three 
points X3 = ±1,0, the value of 3(2:3,71,72) at a generic 
point can roughly be estimated by interpolation, keeping 
the convexity of 5(2:3, 71, 72) with respect to X3 in mind. 
This simple picture is very useful in understanding the 
dependence of 3(2:3,71,72) on various parameters, and 
why the W state is the maximally entangled state with 
respect to GM, as we shall see shortly. 

It is interesting to note that, when X3 > £3 or 

£3 < x 3 3 \ 3(2:3,71,72) is a linear function of 2:3, with 
positive and negative derivatives, respectively. Hence, for 
given 71,72, the minimum of 3(2:3,71,72) is obtained in 

the interval x^ < X3 < x 3 . When 72 = 0, 3(2:3,71,72) 
is an even function of X3 ; its minimum for given 71 is ob- 
tained at £3 = and is equal to ~ . Otherwise, the partial 
derivative of 3(2:3,71, 72) with respect to X3 is positive at 
2:3 = 0; the minimum of 3(2:3,71,72) for given 71,72 is 
obtained in the interval x 3 < X3 < and is smaller than 
\. Setting the derivative of 3(2:3,71,72) with respect to 
X3 to zero leads to a fourth-order polynomial equation 



We have provided analytical methods for deriving 
the GM of symmetric pure multiqubit states with non- 
negative amplitudes in the Dicke basis and that of sym- 
metric pure three-qubit states. Also, we have introduced 
a systematic method for studying the GM of pure three- 
qubit states in virtue of a canonical form of their bipartite 
reduced states. In particular, we have derived explicit an- 
alytical formulae of GM for the family of pure three-qubit 
states one of whose rank-two two-qubit reduced states is a 
convex combination of the maximally entangled state and 
its orthogonal pure state within the rank-two subspace. 
Based on these results, we further proved that the W 
state is the maximally entangled three-qubit state with 
respect to GM. Our studies can simplify the calculation 
of GM and provide a better understanding of multipartite 
entanglement, especially the entanglement in three-qubit 
states. Our results also facilitate the comparison of GM 
with other entanglement measures, like relative entropy 
of entanglement [29;, 46]. Moreover, they may help inves- 
tigate the physical phenomena in multipartite entangled 
systems emerging in condensed matter physics. 
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Appendix A: Derivation of Eqs. (|10H12[I 



When x 1 = x 2 = 0, according to Eqs. ([28]) and ([29]). 

f(a,b,c) = 1 + CCOS71 COS72 + cx 3 sin 71 sin 72 + \w\, 

a(cos72 sin 71 — X3 cos 71 sin 72) 
b{— 23 cos 72 sin 71 + cos 71 sin 72) 
CX3 + X3 cos 71 cos 72 + sin 71 sin 72 



(36) 



We will use the technique in Refs. [36|, |37j to simplify 
the problem. According to the definition in Eq. (JTJ) , 

G 2 (|$)) = max Tr[|$)<$|(|a}<a| ® ® \c)(c\)] 

\a)\b)\c) 

= maxTr[(Tr c |$)(<f|)(| a )(a| ® (32) 

|a>|6) 

where |a),|6) are normalized qubit states. The second 
equality follows from Theorem 1 of E. Jung et at [III, 
which states that any (n— l)-qudit reduced state uniquely 
determines the GM of the original n-qudit pure state, as 
we have mentioned in the second paragraph of Sec. IIIII 

To reduce Eq. (13"2"j) . we use the Bloch sphere represen- 
tation of qubit 22]: 



1 



(I + s p -a), 



(33) 



where the components of <x are three Pauli matrices and 
s p is the Bloch vector. 

Suppose the states \a), \ b) have Bloch vectors s\, s 2 re- 
spectively. Then Eq. (j32[) gives rise to two sets of equa- 
tions: 

ri + Gs 2 = Xxsx, r 2 + Gs 1 =X 2 S2, 

n = Tr[Tr BC (|$X$|)cr], r 2 = Tr[TrAc(|$)($|)o-], 

(34) 

where Ai,A2 are Lagrange multipliers, and the 3x3 
matrix G has elements Gij = Tr[(Trc\$)($\)(cri ® (Tj)\. 
Since the reduced density operators Trecl^K^I an d 
Tr^cl^X^I are identical, one can show that r\ = r 2 = r 
after some algebra. It follows that si = s 2 = s, 
Ai = A2 = A, and Eq. ([3"4"|l reduces to 



r + Gs = As. 



(35) 



The solutions to Eq. (1331) determine the GM of the state 
|$) in Eq. ©. 

Define s — (sin^cos^jSin^sin^jCOsff) with 9 e [0,tt] 
and <p £ [0, 2?r]; then Eq. {35} reduces to Eqs. (fT0[)-(fil 



Appendix B: Derivation of Eq. (|30|) 

In this Appendix, we derive Eq. (j3"0)) . Recall that the 
relevant parameter range is < 72 < 71 < ?, 72+71 < f , 
and —I < x 3 < 1, see Sec. lIII Al To simplify the following 
discussion, we also assume < 71 < 5 and |x 3 | < 1; but 
it turns out that the final result is applicable without this 
restriction. 



According to Sec. IIIIB1 to compute 9(^3,71,72), or 
equivalently, the maximum of /(o,6, c) over the unit 
sphere, we need only to maximize f 2 (c) = f{\l — c 2 , 0, c) 
over the single variable c for — 1 < c < 1. Here f 2 (c) can 
be expressed as follows, 



/2(c) = I + uqc + \Ju\c 2 + 2u 2 c + U3, 
uq = cos 71 cos 72 + X3 sin 71 sin 72 > 0, 
ui = x\ — (COS72 shi7i — X3 COS71 sin72) 2 , 
u 2 = (x 3 cos 71 cos 72 + sin 71 sin 72)2:3, 
u 3 = (sin 7 i) 2 + a;^(cos7i) 2 > 0. (37) 



The four coefficients uq,ui,u 2 ,u 3 in Eq. (|3T[) satisfy the 
following relations, 

u 2 > 1*1, u\ — uiu 3 > 0, u 2 — m > 0; (38) 

these relations are useful in the following discussion. 

To determine the maximum of f 2 (c) for — 1 < c < 1, we 
shall differentiate three cases according to the sign of u\. 
Note that u\ is a quadratic function of x 3 with a positive 
quadratic coefficient, and that it has the following two 
zeros: 



r (i,2) _ cos 72 sin 71 

± I + cos 71 sin 72 ' 



(39) 



which satisfy the inequalities: —I < x^ < < x 1 ^ 1 ' < I 



,.(21 



is equal to —I only if 71 



Case 1: x 3 

uq,u 2 > 0, 



(1) 

x 3 or x 3 



x 



72 = 

(2) 
3 ■ 



1. 

In this case U\ = 0, 



/2(c) = I + u c + \J2u 2 c + 113, 



(40) 



so the maximum of f 2 (c) can only be obtained at c = I . 

Case 2: X3 < x 3 or X3 > x 3 . In this case, 
uo,ui,u 2 > 0, the discriminant of the quadratic func- 
tion u\c 2 + 2M2C + U3 about c is 4(^2 — U1U3) > 0, so the 
quadratic function has two zeros with mean —u 2 /u\ < 0. 
Since the quadratic function must be nonnegative in the 
interval [—1, 1] by definition, both zeros must be smaller 
than or equal to —I. In the interval [—1,1], this quadratic 
function and /2(c) are both strictly increasing, so the 
maximum of f 2 (c) can only be obtained at c = I. 

Case 3: x^ < X3 < x^\ In this case u\ < 0, the 
quadratic function u\C 2 + 2u 2 c + 1*3 is positive between 
its two zeros. One zero is smaller than or equal to — 1, 
and the other larger than or equal to I. To determine 
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the maximum of /2(c), we take the first and the second 
derivatives of /2(c): 



UQ 



u\c + u 2 



y/uiC 2 + 2u 2 c + u 3 



U1U3 — U- 2 



(u^c 2 +2u 2 c + u 3 ) 3 /' 



<0, 



(41) 



where the last inequality follows from Eq. ([55]). There is 
only one solution to the equation /2(c) = 0, 

c = — \x 3 — (cos 72 sin 71 — x 3 cos 71 sin 72) 2 ] 
x [^3(^3 cos 71 cos 72 + sin 71 sin 72) 
+ sin 71 (cos 71 cos 72 + £3 sin 71 sin 72) 
x (sin 71 — x 3 COS71 tan 72)] > 0. (42) 

Since the second derivative of /2(c) is always negative, c is 
the global maximum of the function f 2 (c) in the interval 
where it is real valued. Restricted to the interval [—1,1], 
the maximum of /2(c) is obtained at c if c < 1 and at 
c = 1 otherwise. In both cases, the maximum point is 
unique. Hence, it remains to determine when c > 1 and 
when c < 1. 

( 1 2) 

After some algebra, one can show that x 3 ' defined 

in Eq. and zf' 4 ' defined in Eq. (|3Tj! satisfy the fol- 
lowing inequalities, 



1< x 



(2) 



< Xg 3 '' < < £3*' < X3 



,.(4) 



(1) 



If 7i + 72 = § , then 1E3 
following relation, 



(3) 



,.(' 2 ) 



< 1. (43) 
1, and c satisfies the 



where in deriving the last equality, we have noticed that 

{(3) 
X3 + %3 cos 71 cos 72 + sin 71 sin 72 < 0, x 3 < x 3 , 
x 3 + x 3 cos 71 cos 72 + sin ji sin 72 > 0, 2; 3 > x 3 4 ^ . 

Now Eq. ([3"U| follows immediately when < 71 < ? and 
|x 3 | < 1; recall that g(x 3 , 71,72) = 3 max_i< c <i / 2 (c). 
It is straightforward to verify that the formula is also 
valid in the special cases \x 3 \ = 1 or 71 = 0, 77, hence the 
derivation is complete. 



Appendix C: Proof of Theorem [Tl 



To prove Theorem [T] in Sec. IIII CI it suffices to show 
that the global minimum of g(p r k2) is obtained at the 
two-qubit reduced state of the W state. 

In the relevant parameter range < 72 < 71 < \, 
71+72 < § , by taking its derivative with respect to 72 in 
Eq. (|30|) . one can show that, for given 71, x 3 , g(x 3 , 71, 72) 
is monotonically decreasing with 72 for x 3 < and mono- 
tonically increasing with 72 for x 3 > 0. Assuming x 3 < 0, 
where the global minimum point should satisfy according 
to Sec. IIII Cl then 3(^3,71,72) is monotonically decreas- 
ing with 72. Hence, either 71 = 72 or 71 + 72 = \ at the 
global minimum of g(x 3 , 71, 72). We shall show that the 
unique minimum is obtained at the two-qubit reduced 
state of the W state in both cases. 



1. special case: 71 + 72 = f 




i< x 3 < 4 4) , 



..(4) 



< x 3 < x 



(1) 



(44) 



If 71 + 72 < § , then — 1 < x^ < x^' , and c satisfies the 
following relation, 



c< 1 
c> 1 



4 <x 3 < 4 4) , 

(2) ^ (3) 

x 3 < x 3 < x 3 



(4) . ^ (i) (45) 
or x\ < x 3 < x\ . 



According to the observations in the above three cases, 

(3) (4 

\ or x\ 



if — 1 < x 3 < Xg 3 ' or x 3 4 ^ < x 3 < 1, the maximum of /2(c) 



is obtained at 1; if — x 3 °' < x 3 < x 3 4 \ the maximum is 
obtained at c. The maximum point is unique in both 
cases. The values of /2(c) at 1 and c are respectively 
given by 

/»(*) = 

2 ( 1 — x 3 ) sin 71 cos 72 (cos 72 sin 71 — x 3 cos 71 sin 72 ) 

— \x 3 — (cos 72 sin 71 — x 3 cos 71 sin 72 ) 2 ] 
/ 2 (1) = 1 + cos 71 cos 72 + x 3 sin 71 sin 72 

\x 3 + x 3 cos 71 cos 72 + sin 71 sin 72 1 

(1-X 3 )[1+ C0S( 7 1 + 72)] , X 3 < , ( 4 g) 

(1 +x 3 )[l +cos(7i -72)], a;3>4 4) , 



If 72 = f - 7i (f < 7i < f ), Prk2 is supported on the 
symmetrical subspace, according to Eqs. (|22p -(|24 |) . In 
this case, Eq. ([30]) reduces to 



{1 (l+Z3)a;3(cos7i) 2 
2 -l+3K3 + (l+x 3 ) cos(2 7 i) : 
(l+x 3 )[l+ s in(2 7l )] 
4 ' 



T<£3<4 4) : 



<x 3 <l, 



(47) 



where 



l-V2sin(2 7 i + f) 
3 + V2sin(2 7 i + f)' 



(48) 



(7(0:3, 71, 5 — 71) is equal to ^ &t x 3 — 0, —1, independent 
of 71; <7(x3,7i, f — 71) is monotonically increasing with 
71 for —1 < x 3 < 0, and monotonically decreasing for 
< x 3 < 1 (see also Fig. p. 

To determine the maximum of g(x 3 , 71, ^ — 71) for 
given 71, we can set its derivative with respect to x 3 
to (in the interval — 1 < x 3 < x^), which leads to the 
following quadratic equation over x 3 , 



[3 + cos(2 7 i)]4 + [-2- 



2cos(27i)]x 3 +cos(27i) = 1. 

(49) 
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Given f < 71 < §, the only solution with modulus less 
than or equal to 1 is 



where 



(5) _ 2 sin 71 (sin 71 - \/2) 
Cs ~ 3 + cos(2 7 i) 



(50) 



,(3) 



-(tan7i) 2 . 



(54) 



The minimum of 3(2:3, 71, § — 71) for given 71 is 



I (5) 
3(4 : 



71 



7T 

2 



71 



[l + cos(27i) + -\/2sin7i] 5 



[3 + cos(2 7 i)] 



(51) 



One can show that g(x^\ 71, £ — 71) is monotonically 
increasing with respect to 71 by taking its derivative with 
respect to 71; hence, its minimum is obtained at 71 = j. 

(5) 1 

At this minimum point, 71 = 72 = tj 2:3 = £3 = — -3=, 
and <?(—§) J, f) — |. This minimum is also the global 
minimum of <?(p r k2)- 

The rank-two state corresponding to this minimum is 
exactly the two-qubit reduced state of the W state, more- 
over, up to local unitary transformations, the W state is 
the only pure three-qubit state with this rank-two state 
as a two-qubit reduced state. To see this, recall that the 
two-qubit reduced state of I W) = ^(|100)+|010)+|001)) 



p rk2 (W) = I|00)(00| + ^ + )(^+| (52) 



with \ip^ 



1 



(|01) + 1 10) ) , which is a convex combina- 



tion of a pure product state and a orthogonal Bell state, 
thus 71 = 72 = f > according to Ref. [45]. In this rank- 
two subspace, the Bloch vectors of the two states |00) 
and \ip + ) are (0,0,1) and (0,0,-1) respectively, and the 
Bloch vector of the state p r k2(W) is exactly (0,0, — |). 



It is interesting to note that 17(23, 71, 71) is independent of 
71 when < X3 < 1. If 2:3 < 0, 3(2:3, 71, 71) is monotoni- 
cally decreasing with 71 . Hence, 71 = j at its minimum, 
the corresponding rank-two subspace is then symmetric. 
According to the result on symmetric states in the pre- 
vious subsection, the unique minimum of 3(2:3,71,71) is 
also obtained at 71 = 72 = 2:3 = — 5. 

We have shown that the unique minimum of 
3(2:3,71,72) is obtained at 71=72 = |, X3 = -±, 
and that the corresponding state is the two-qubit re- 
duced state of the W state. This minimum is also the 
global minimum of g(p r \a)- To prove that the min- 
imum is unique among all two-qubit rank-two states, 
it remains to show that it is unique in the rank-two 
subspace with 71 = 72 = ?■ It suffices to ver- 
ify that 3(2:1,2:2, 3 , 4 , 

x\ > (here we write 3(2:1, x 2 , 2:3, 71, 72) for 3(^2))- 
Due to the rotational symmetry of 3(2:1, X2, 2:3, 71, 72) 
about the 2:3 axis discussed in Sec. IIII A[ this is true if 
3(2;i,0,-i,f,f) > 3(0,0, -i,f,f) for x x > 0. Accord- 
ing to Eq. ([35]), when 71 = 72 = f , 2:2 = 0, 2:3 = -i, 



J > ff (0,0,-|,f,f) for x\ + 



■ [5 + 32:1 + v/9 + 32;i(10 + 9xi)" 



(55) 



2. special case: 72 = 71 < \ 



hence, 



In this case, Eq. ([30)1 reduces to 

3(^3,71,71) = 

^(cos 7 i) 2 , 

-(l-a: 3 ) 2 (l+a: 3 )[ S in(27 1 )] 2 
-l+x a {2+7x 3 ) + {l-x 3 y 2 cos(4 7 i) 
(l+x 3 ) 



i<x 3 <4 3) , 



,(3) 



< 2:3 < 0, 



< x 3 < 1, 



3(^1,0,- 



1 ,{2y/2 1 
'3 



1 _ -)>-/( 

3' 4' 4/ ~ A J \ 3 

4 / 1 7T 7T 

> 9 =3(0,0, 



3' 4' 4 



(53) This completes the proof of Theorem [T] 



(56) 
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